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ABSTRACT 

This  paper  is  devoted  to  the  presentation  of  the  general  framework  and  the  initial  results 
of  a  joint  effort  to  derive  novel  research  tools  and  easy  to  use  software  to  analyze  and  model 
turbulence  and  transition. 

After  a  brief  review  of  the  issues  and  a  summary  of  some  basic  properties  of  wavelets,  we 
present  our  preliminary  results.  Both  the  technical  aspects  of  the  implementation  and  the 
physical  conclusions  reached  at  this  time  are  discussed. 

Current  developments  are  summarized  in  the  last  section. 


'Research  was  supported  in  part  by  the  National  Aeronautics  and  Space  Administration  under  NASA 
Contract  No.  NASI-18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665 
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1.  INTRODUCTION 


The  goal  of  this  work  is,  through  a  common  effort  of  people  from  various  fields  (Turbulence 
theory,  Experimental  and  Numerical  simulation,  Signal  processing,  Numerical  analysis,  Har¬ 
monic  analysis),  to  generate  innovative  research  tools  devoted  to  the  analysis,  understanding 
and  modeling  of  turbulence. 

The  starting  point  of  this  program  is  the  wavelet  decomposition  and  the  initial  algorithms 
derived  by  some  of  the  authors  (Liandrat  and  Moret-Bailly  1990,  Liandrat  and  Tchamitchian 
1990,  Moret-Bailly  et  al.  1991)  to  study  and  model  turbulence,  and  more  generally,  nonlinear 
phenomena.  Compared  with  more  classical  analysis  tools,  e.g.  Fourier  Transforms,  this  new 
family  of  algorithms  is  very  flexible.  However,  it  rapidly  became  apparent  that  before  these 
algorithms  become  really  useful  to  the  fluid  dynamics  scientific  community,  they  must  first 
be  organized,  and  translated  into  an  understandable  and  easily  usable  form. 

In  this  interim  report,  we  present  the  first  kernel  of  a  program  devoted  to  this  task.  The 
first  applications  and  the  ultimate  goal  of  this  work  are  also  briefly  presented.  We  wish  to 
emphasize  that  the  goal  of  the  work  is  to  greatly  extend  the  now  classical  wavelet  decom¬ 
position  algorithms.  A  hierarchy  of  tools  is  being  developed  that  will  provide  quantitative 
results  together  with  the  basic  elements  for  the  modeling  of  variables  within  the  context  of 
nonlinear  dynamics,  transition  and  turbulence.  Currently,  several  algorithms  developed  for 
the  study  of  boundary-layer  transition  on  a  rotating  disk  have  been  implemented  and  are 
being  tested  in  other  configurations. 

Continuous  interaction  between  all  the  collaborators  of  this  program  (i.e.  the  authors 
of  this  report)  should  not  only  bring  the  already  existing  algorithms  to  a  level  of  easy 
applicability,  but  will  also  provide  two  additional  benefits.  First,  a  subset  of  these  algorithms 
will  evolve  into  methodologies  which  will  be  directly  related  to  the  modeling  of  the  studied 
phenomena.  An  example  of  this  is  the  relation  between  wavelet  analysis  tools  and  Large 
Eddy  Simulation.  The  second  benefit  of  this  interaction  will  be  the  give  birth  to  a  new 
generation  of  algorithms,  more  powerful  than  the  last. 

2.  REVIEW  OF  WAVELET  TRANSFORM  BASICS  AND  NOTATIONS 

Since  the  basic  motivation  of  this  work  is  the  wavelet  decomposition,  some  basic  proper¬ 
ties  of  one-dimensional  wavelets,  together  with  the  conventions  used  in  the  report  are  first 
reviewed. 

A  complete  review  on  the  current  state  of  wavelet  theory  is  not  available,  but  various 
tutorials,  book  conference  or  reviews  can  be  found  in  the  literature  (Combes  et  al.  1989, 
Rioul  and  Vetterli  1991,  Farge  1992) 


From  an  analysis  function  xj)(x)  (Figure  1),  the  family  of  functions  used  to  decompose  a 
signal  u( x)  is  generated  by  dilations  and  translations  following  the  formula: 

(1)  ^ab{x)  =  ~  l/>(- - -) 

y/a  a 

where  («,  b )  belongs  to  a  2D  continuous  or  discrete  space,  (a,  b )  belongs  to  R*+  x  R  in  the 
so-called  continuous  case  and  to  a  countable  subspace  in  the  discrete  case.  If  in  addition, 
xpab  form  an  orthogonal  basis,  one  has:  a  =  2J,  b  =  k2~^  with  ( j,  k)  £  2Z  2  and  the  dilated 
and  translated  wavelet  is  denoted  x bjk  instead  of  xj)ab. 

One  defines  the  wavelet  coefficients  ( T(b,a )  or  7)*)  of  the  function  u(x )  as: 


or  equivalently  as: 


/+oo 

u(x)xp*ab(x)dx 

-OO 

/+oo 

u(v)fcb(oj)dw 

-oo 


where  u’(x)  and  u(u;)  stand  respectively  for  the  conjugate  and  the  Fourier  transform  of  the 
function  u(z). 

One  recovers  the  original  function  u(x)  through  summation  formulas  coming  from  the 
identity  decomposition  (see  Combes  et  al.  (1989)  for  details). 

In  the  continuous  version,  one  recovers  the  original  function  from 


1  f+oo  r+oo  dac 

/  /  T(b,a)xbah(x )  — 

Jo  j- oo  o 


C(rP) 


and  its  energy  ( L 2  norm)  from 


/+oo  1  r-  OO  r+OO  (Jnr 

u(x)u‘(x)dx  =  /  /  T(b,  a.)T(b,  a)‘  - 

-oo  Clxb)  Jo  J- oo  a* 
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where  C(xl>)  is  a  constant  depending  only  on  xji(x). 

For  the  orthogonal  decomposition  one  gets: 

(6)  u(x)  = 

jk 


/+oo  _ _ 

u(x)u(x)dx  =  Y,Tv 

-OO  •  . 


The  wavelet  theory  is  ’’constructive”  in  the  sense  that  the  construction  of  the  wavelets 
brings  along  with  it  the  algorithms  needed  to  compute  t  he  wavelet  coefficients  of  the  decom¬ 
position.  Most  of  these  algorithms  have  an  asymptotic  operation  count  at  most  equal  to  that 


of  the  Feist  Fourier  Transform  which  is  0(N  log  N).  If  the  wavelet  is  of  compact  support, 
O(N)  algorithms  exist.  Here,  N  stands  for  the  number  of  points  in  the  signal).  Details  on 
the  algorithms  can  be  found  in  the  book  of  Combes  et  al.  (1987)  and  in  the  papers  of  S. 
Mallat  (1988)  and  Perrier  and  Basdevant  (1989). 

As  stated  in  the  introduction,  the  wavelet  plane  (i.e  the  family  of  the  wavelet  coefficients 
Tab)  is  assumed  to  be  the  zero  level  output  of  the  algorithms.  This  2D  field  becomes  the 
starting  point  for  the  first  level  algorithms. 

3.  SPECIFIC  PROBLEMS  OF  TURBULENCE  AMENABLE  TO  WAVELET 
ANALYSIS 

If  one  had  to  summarize  the  initial  motivation  for  the  use  of  wavelets  in  the  field  of 
turbulence,  localization  and  scaling  would  be  two  key  words  for  both  wavelets  and  turbulence. 
These  two  properties  can  be  used  to  good  advantage  to  help  decompose  turbulence  into  its 
natural  scales  without  loss  of  spatial  information,  and  also  to  compress  turbulent  data  by 
throwing  away  information  which  is  not  directly  relevant  to  the  observed  turbulent  signals. 

But,  as  stated  earlier,  the  fundamental  question  addressed  here  is  the  following: 

What  quantitative  information  can  be  extracted  from  the 
wavelet  decomposition  of  a  data  field ? 

Within  that  context,  a  good  choice  of  the  wavelet  family  is  important  since  the  ’’quanti¬ 
tative  information”  must  refer  to  the  analyzed  signal  and  not  to  the  analyzing  wavelet.  In 
that  regard,  it  is  known  that  localization  in  both  physical  and  Fourier  spaces  is  necessary. 
Numerically,  spline  wavelets  of  sufficiently  high  degree  (say  m  >  4)  are  satisfactory  (Liandrat 
and  Tchamitchian  1990). 

After  a  choice  for  the  family  of  wavelets  is  made,  we  turn  our  attention  to  the  wavelet 
plane.  Energetic  information  (Moret  Bailly  et  al.  1991)  is  first  derived  and  a  new  plane 
of  coefficients  is  built.  This  energy  plane  now  forms  a  new  platform  upon  which  second 
level  algorithms  are  constructed.  Different  representations  and  decompositions  of  this  plane 
are  then  available.  Depending  on  there  definitions,  they  can  refer  for  instance  to  coherent 
structures  or  optimal  decomposition  of  the  signal.  These  issues  are  at  the  cutting  edge  of 
data  interpretation  in  turbulence. 

Finally,  it  appears  that  a  large  number  of  specific  problems  should  be  addressed  by  this 
program.  Among  them  are: 
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1.  The  study  of  the  relations  between  the  global  spectrum  (or  scale  decomposition)  and 
singularity  distributions  in  a  turbulent  signal.  This  problem  relates  to  the  fractal  or 
multifractal  properties  of  turbulence  as  well  as  the  notions  of  complex  singularities  in 
the  theory  of  partial  differential  equations  (Sulem  et  al.  1983). 

2.  The  characterization  of  processes  undergoing  transition  using  local  scale  decomposi¬ 
tions  and  as  a  consequence,  of  the  development  of  models  for  turbulence. 

3.  The  definition  and  extraction  of  typical  (coherent)  events  and  the  quantification  of  the 
vortex  modeling  of  turbulence. 

Answers  to  these  issues  would  lead  to  a  better  understanding  of  the  mechanisms  intrinsic 
to  the  dynamics  of  the  flow.  This  understanding  will  in  turn  lead  to  a  new  generation  of 
modeling  concepts  for  the  turbulence  cascade,  and  for  transition  and  turbulence  in  general. 

4.  ALGORITHMS  ALREADY  IMPLEMENTED 

In  this  section,  the  algorithms  that  have  already  been  implemented  along  with  some 
significant  results  are  presented.  One  should  refer  to  the  quoted  papers  for  more  details. 

4.1.  Basic  wavelet  decomposition  algorithms  (Zero  Level) 

Two  kinds  of  decomposition  algorithms  have  been  implemented  that  provide: 

-  A  decomposition  on  periodic,  even  order  spline  orthogonal  wavelets  (for  N  =  2P  data 
points,  a  =  2~Jand  b  =  k2~J  where  0  <  j  <  p  —  1  and  0  <  k  <  2J  —  1) 

-  A  decomposition  on  periodic,  even  order  spline  wavelets  (for  N  =  2P  data  points, 

a  =  =  k  where  1  <  n  <  y  and  0  <  k  <  N  —  1)  (see  [Perrier  and  Basdevant  1989]) 

This  last  algorithm  provides  a  fully  redundant  decomposition  in  the  position  parameter 
b.  The  ratio  of  the  number  of  voices  per  octave  can  be  adapted  using  available  interpolation 
routines. 

In  both  algorithms,  the  scale  limits  can  be  arbitrary  chosen. 


4.2  First  Level  Algorithms 

Starting  from  the  wavelet  plane,  the  energy  plane  is  constructed.  Following  equation  (5), 
the  energy  at  the  point  (b,  a)  is  (Moret-Bailly  et  al.  1991): 


(8) 


E(b,a)  = 


1 


|7\M|2 


C(0) 


r\m 

J  a 


1 

-x 

a 


'b'  -b'' 


dti 


where  \(x)  is  a  bump  function  that  represents  the  envelope  of  r['(x)  and  satisfies 
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J  X{x)dx= 1.  This  leads  to  J  E(b,a)dadb  =  ||u||2. 

Using  a  different  normalization  at  each  point,  the  local  energy  density  probability  (also 
called  scale  decomposition)  is  defined  as 


(9) 


Db(a) 


E(b,  a) 

J  E(b,a)da 


4.3  Second  Level  Algorithms 

The  analysis  of  the  energy  plane  E(b,a)  in  terms  of  the  characterization  of  local  or 
averaged  scale  decomposition.  Over  a  subspace  S  of  points,  the  averaged  scale  decomposition 
is  given  by 


Ds(a)  = 


Eb6sg(M) 

/ (H  E{b,a))da 

J  L  C 


Starting  from  this  scale  decomposition,  one  defines  a  mean  scale  as: 


as  =  J  log (a)Ds(a)da 


where  the  subscript  s  referes  to  a  point  b  or  a  subset  S.  Notice  that  the  integration  is  made 
in  function  of  log(a)  and  not  in  function  of  a.  Indeed,  log(a)  is  the  scale  that  corresponds 
to  the  dilation  factor  a. 

A  normalized  standard  deviation  w  is  then  introduced  as: 


ws  =  J{log(a)  -  as)2 


The  mean  scale  a3  is  the  scale  around  which  the  active  scales  gather  at  point  b  (or  in 
the  subspace  S )  while  the  normalized  standard  deviation  ws  quantifies  the  dispersion  of  the 
active  scales  around  as. 

The  scale  decomposition  averaged  over  the  entire  energy  plane  is  easily  related  to  the 
Fourier  energy  spectrum  thanks  to  equation  (3).  A  case  of  special  interest  for  turbulence 
occurs  when  the  energy  spectrum  exhibits  a  u>2a  behaviour  at  large  frequency.  Then,  the 
averaged  scale  decomposition  behaves  as  a_2o_2.  (see  example  5.2). 

The  local  scale  decomposition  characterizes  the  scaling  behaviour  of  a  function  at  a  given 
point  (Tchamitchian  and  Holschmeider  1989).  For  a  Holder  exponent  a  at  point  b  the  scale 
decomposition  at  the  point  b  behaves  as  a2o+1.  (see  example  5.3). 


4.4  Ultimate  Level  Algorithms 


Derived  for  a  specific  application,  these  algorithms  make  the  connection  with  the  physical 
modeling  of  the  phenomenon  under  study.  The  existing  algorithms  are  being  applied  to  the 
study  of  the  transition  on  a  rotating  disk  and  to  the  analysis  and  characterization  of  wall 
structures.  One  must  refer  to  Moret-Bailly  et  al.  (1990)  for  details. 

5.  FIRST  RESULTS 

In  the  following  section  the  graphical  environment  of  the  existing  version  of  the  code  is 
presented  along  with  some  illustrations  of  the  algorithm  outputs  at  each  level. 

5.1  Test  case  of  modulated  waves  (see  Figure  2) 

For  this  test  case  (Figure  2. a),  the  local  mean  scale  (Figure  2.b),  standard  deviation 
(Figure  2.c)  and  local  scale  decomposition  (Figure  2.d)  plots  are  presented.  The  two  flat 
zones  of  the  mean  scale  decomposition  represent  the  local  waves,  while  the  local  peaks  in 
the  standard  deviation  reveals  the  singular  regions.  The  different  local  scale  decompositions 
can  be  used  to  characterize  each  part  of  the  signal. 

5.2  ui2a  energy  spectrum  signal 

A  signal  computed  from  a  random  phase  u j~*  energy  spectrum  has  been  constructed 
and  wavelet  transformed  (Figure  3. a).  The  averaged  scale  decomposition  (Figure  3.c)  is 
compared  to  the  Fourier  spectrum  (Figure  3  b).  The  local  scale  decompositions  (Figure  3.d) 
clearly  show  that  no  local  scaling  occurs  in  the  signal  (see  Figure  3. a). 

5.3  Test  case  of  local  singularities 

The  analyzed  signal  exhibits  a  x2  local  behaviour  around  i  =  512  (Figure  4. a).  This 
local  scaling  is  revealed  by  the  local  scale  decompositions  (Figure  4.b)  which  gives  an  a 5 
spectrum. 

5.4  Rotating  disk  boundary  layer  velocity  signal 

Figures  5  show  some  results  obtained  from  the  analysis  of  a  hot  film  velocimeter  signal 
(Figure  5. a)  in  a  transitional  boundary  layer  over  a  rotating  disk.  As  a  first  step,  the 
variations  of  the  standard  deviation  (Figures  5.b  and  5.c)  are  used  to  discriminate  turbulent 
bursts  from  smoothly  oscillating  parts  of  the  signal.  Then,  a  full  decomposition  of  the  signal 
parametrized  by  the  value  of  the  standard  deviation  is  performed. 

A  comparison  of  this  decomposition  at  different  Reynolds  numbers  is  used  to  modelize 
the  transition  from  an  oscillating  regime  all  the  way  to  turbulence  (Moret-Bailly  et  al.  1991 ). 
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5.5  Analytical  signals  to  model  the  compressible  boundary  layer 

A  computer  code  that  generates  transitional  signals  has  been  written  for  use  in  connec¬ 
tion  with  a  digital  signal  processing  (DSP)  system  under  development  at  NASA  Langley. 
This  system  will  be  part  of  hypersonic  flight  tests  (Bertelrud  1991).  The  output  signal  is  de¬ 
terministic  and  has  built  into  it  intermittency,  simulated  anomalies  and  a  variety  of  periodic 
or  turbulent  quantities  undergoing  a  series  of  discrete  oscillations.  The  signal  chosen  for  the 
current  analysis  consisted  of  either 

1.  A  single  frequency  in  the  laminar  regime  (like  a  wind  tunnel  fan  might  yield) 

2.  A  series  of  frequencies  with  an  amplitude  distribution  corresponding  to  the  —1/3  and 
—5/3  slopes. 

The  full  signal  contains  a  DC  as  well  as  an  AC  component.  However,  the  DC  component 
was  removed  from  the  signal  displayed  in  Figure  6. a. 

The  analysis  is  performed  in  term  of  mean  scale  (Figure  6.b)  and  standard  deviation 
(Figure  6.c). 

The  intermittency  function  that  discriminates  the  different  parts  of  the  signal  (6.d)  is 
obtained  from  the  study  of  the  standard  deviation.  A  comparison  between  the  averaged 
scale  decompositions  in  each  part  of  the  signal  with  an  adapted  windowed  Fourier  transform 
is  presented  in  Figure  (6.e).  Here,  the  averaged  scale  decomposition  is  plotted  versus  the 
frequency  instead  of  the  dilation  parameter  a  (one  reminds  that  the  bump  shape  of  the 
Fourier  transfer,  i  of  a  wavelet  allows  to  associate  to  each  value  of  a  a  frequency  oo). 

6.  TECHNICAL  IMPLEMENTATION 

6.1  Goals 

The  technical  implementation  of  the  algorithms  and  the  associated  graphics  should  be 
clear,  interactive,  well-suited  to  incorporate  new  algorithms,  and  easily  portable  across  com¬ 
puter  platforms.  Among  the  possible  solutions  to  this  challenge,  three  types  of  interfaces  are 
under  study: 

-  A  full  graphic  user  interface  with  buttons,  sliders,  ...  . 

-  An  interconnected  family  of  modules 

-A  graphic  system  governed  by  a  specific  command  language. 

6.2  Existing  version 

At  this  time,  our  prototype  has  the  first  type  of  user  interface.  Typical  screens  can  be 
seen  on  figure  7  with  the  different  windows  standing  for  the  signal  representation  (1),  the 
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wavelet  or  energy  plane  (2),  the  local  or  mean  scale  decomposition  (3),  the  second  level 
output  window  (4)  and  the  choice  and  control  window  (5).  Although  this  software  is  already 
quite  satisfactory,  it  does  not  satisfy  all  the  previously  stated  objectives.  However,  we  think 
that  this  version  can  already  help  answer  basic  questions  and  will  serve  as  an  elementary 
module  for  the  next  release  of  the  program. 

7.  FUTURI  DEVELOPMENTS  AND  CONCLUSIONS 

Together  with  the  technical  implemeiA ation,  new  algorithms  have  been  constructed  and 
tested.  On  the  one  hand,  the  second  level  algorithms  described  in  this  report  are  being  refined 
by  the  introduction  of  a  better  description  of  the  local  or  averaged  scale  decompositions.  An 
algorithm  devoted  to  the  detection  and  characterization  of  complex  singularities  is  being 
implemented  and  tested. 

Third  level  algorithms  based  on  the  optimal  L2  decomposition  (Karhunen-Loeve  decom¬ 
position)  of  the  energy  plane  are  under  study  in  the  framework  of  coherent  structures  and 
transition  modeling. 

On  the  other  hand,  among  the  algorithms  already  available,  the  class  of  ridge  and  skeleton 
algorithms  derived  by  Toressani  and  Tchamitchian  (1991)  and  the  class  of  fractal  character¬ 
ization  algorithms  developed  by  Arneodo  et  al.  (1991)  should  be  implemented  numerically 
in  the  near  future.  They  clearly  have  applications  (even  if  it  is  not  yet  completely  demon¬ 
strated)  in  the  fields  of  turbulence  and  transition. 
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Figures  6:  Model  of  compressible  boundary  layer  signal 
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Local  scale  decomposition  Areraged  scale  decomposition 


Figure  7:  Screen  from  graphic  user  interface. 
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